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00 ■ Abstract 

fSJ ■ We study the impact of minimal non-supersymmetric models of resonant leptogenesis on charged 
' lepton flavour violation and the neutrino mixing angle ^13. Possible low-scale flavour realisations 
of resonant r-, /i- and e- leptogenesis provide very distinct and predictive frameworks to explain 
the observed baryon asymmetry in the Universe by sphaleron conversion of an individual r-, /x- 
and e-lepton-number asymmetry which gets resonantly enhanced via out-of-equilibrium decays 
of nearly degenerate heavy Major ana neutrinos. Based on approximate flavour symmetries, we 
construct viable scenarios of resonant r-, /x- and e- leptogenesis compatible with universal right- 
handed neutrino masses at the GUT scale, where the required heavy-neutrino mass splittings 
are generated radiatively. The heavy Majorana neutrinos in such scenarios can be as light as 
100 GeV and their couplings to two of the charged leptons may be large. In particular, we 
explicitly demonstrate the compelling role that the three heavy Majorana neutrinos play, in order 
to obtain successful leptogenesis and experimentally testable rates for lepton flavour violating 
processes, such as /i — t- 67 and // — )• e conversion in nuclei. 
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1 Introduction 



The observed baryon asymmetry in the Universe (BAU), which amounts to a baryon-to- 
photon ratio of number densities rjB ~ 6.2 x 10~^° [lll2], provides one of the strongest pieces 
of evidence for physics beyond the Standard Model (SM) One interesting scenario for 
explaining the BAU is leptogenesis [4j . Leptogenesis does have a profound link to neutrinos 
and the origin of their extraordinary small masses. In particular, the famous seesaw mech- 
anism |5j can give rise to their small observed masses through the presence of superheavy 
Majorana neutrinos close to the Grand Unified Theory (GUT) scale, Mqut ~ 2x 10^^ GeV. 
These GUT-scale heavy neutrinos, being singlets under the SM gauge group, may possess 
large Majorana masses that violate lepton number (L) by two units. In an expanding 
Friedmann-Lemaitre-Robertson-Walker (FLRW) Universe, the heavy Majorana neutrinos 
can decay out of equilibrium and produce a net leptonic asymmetry. The so-produced 
leptonic asymmetry gets rapidly reprocessed into the observed BAU [1] by equilibrated 
{B + L)-violating sphaleron interactions [5] . 

A potentially interesting alternative to GUT-scale leptogenesis is the framework of 
low-scale resonant leptogenesis (RL) [TJE] . Within this framework, the lowering of the scale 
may rely on a dynamical mechanism, in which heavy-neutrino self-energy effects [9] on the 
leptonic asymmetry become dominant [10] and get resonantly enhanced [7] , when a pair of 
heavy Majorana neutrinos has a mass difference comparable to the heavy neutrino decay 
widths. As a consequence of thermal RL, the heavy Majorana mass scale can be as low as 
~ 100 GeV, while maintaining agreement with the solar and atmospheric neutrino data |8] . 
One of the advantages of low-scale RL is that the reheating temperature Treh resulting from 
inflaton decays does not need to be very high, e.g. T^eh ~ 1^10 TeV JTT\, thereby avoiding 
comfortably the overproduction of gravitinos in supersymmetric models whose late decays 
may cause dissociation of the light elements during the nucleosynthesis era [T2] . 

Flavour effects due to heavy-neutrino Yukawa couplings play an important role in 
models of RL [T3l[Tl]. In particular, in [13] a scenario was put forward, called resonant 
r-leptogenesis (RrL), in which the BAU originates from a r- lepton asymmetry, resonantly 
produced by quasi-in-equilibrium decays of heavy Majorana neutrinos. This mechanism 
makes use of the property that sphalerons preserve the individual quantum numbers — 
Le,^i,T [T514T8] . In a RrL model, the generated excess in the L^. number will be converted 
into the observed BAU, provided the L^-violating reactions are not strong enough to wash 
out such an excess. In such a scenario, the heavy Majorana neutrinos can be as light 
as 100 GeV and have sizeable couplings to two of the charged leptons, specifically to 
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electrons and muons. Consequently, depending on the flavour dynamics of heavy neutrino 
Yukawa coupling effects, phenomenologically testable models of RL can be built that could 
be probed at the LHC or in low energy experiments of lepton number violation (LNV) 
and lepton flavour violation (LFV). For instance, observables of particular interest are the 
neutrinoless double-beta {0^(3(3) decay of heavy nuclei [19], the photonic decay /i — )■ 67 
analyzed by the MEG experiment ^20j and the coherent — )■ e conversion in nuclei to be 
looked for in the planned COMET/PRISM experiment [2T] . 

In this paper we study all possible alternatives to RrL, including the minimal models 
of resonant /x-leptogenesis (R/xL) and resonant e-leptogenesis (ReL). Collectively, we refer 
to these three different lepton-flavour realisations of RL as R£L. We assume that the R£L 
models have an S0(3)-symmetric heavy neutrino mass spectrum at the GUT scale, with all 
heavy Majorana neutrinos being exactly degenerate at this scale. We consider a minimal 
non-supersymmetric framework, in which the heavy-neutrino mass splittings required for 
successful RL are generated radiatively [22] and can therefore be naturally comparable to 
the decay widths of the heavy neutrinos. Since all charged-lepton Yukawa couplings are 
in thermal equilibrium at temperatures T < 10 TeV [23], we consider a flavour diagonal 
basis for these couplings, while setting up the Boltzmann equations (BEs). In addition, 
we include the flavour effects due to individual heavy-neutrino Yukawa couplings [T3|[2l]. 
which can have a dramatic impact on the predictions for the BAU in RL models [T3|[Ti]. 

The layout of the paper is as follows. Section |2] describes the basic structure of the 
minimal SM with three heavy Majorana neutrinos and introduces the flavour symmetries 
needed to realize the different lepton-flavour scenarios associated with KiL. As mentioned 
above, we assume a SO (3) symmetric heavy- neutrino sector at the GUT scale and calculate 
the renormalization-group (RG) effects on the mass spectrum of the electroweak-scale heavy 
Majorana neutrinos and their Yukawa couplings to charged leptons. Taking the light- 
neutrino oscillation data into account, we are able to determine most of the theoretical 
parameters of the RiL models. In Section [3l we present analytic results and predictions 
for LFV observables in the three different BIL models. Section H] briefly reviews the basic 
framework of RL and presents the BEs, upon which our numerical analysis is based. We 
also clarify the necessity of having at least three heavy Majorana neutrinos in RL models 
in order to obtain experimentally testable LFV. Section [5] presents numerical estimates of 
representative R£L models and their impact on the neutrino mixing angle ^13. Finally, 
Section El summarizes our conclusions. 



3 



2 Flavour Models of Minimal Resonant Leptogenesis 



In this section, we describe the basic theoretical framework underlying the different flavour 
models of minimal RL. In particular, our interest is in scenarios in which the BAU is 
generated by the production of an individual lepton number [13]. For definiteness, we first 
consider a minimal model for RrL in Section [2TT1 and then generalize to the other two cases 
R/iL and ReL in Sections 12.21 and 12.31 respectively. 

The leptonic Yukawa and Majorana sectors of the SM symmetrically extended with 
one singlet right-handed neutrino um per i family (with i = 1,2,3 = e,fi,T) are given by 
the Lagrangian 

- Cy,m = L^hUn + L^h^UR + u'imMiyR + H.c, (2.1) 

where $ is the SM Higgs doublet and $ = ir2$* its isospin conjugate. Moreover, we 
have suppressed the generation index i from the left-handed doublets Lj = (z/j^,, Ul)'^, 
the right-handed charged leptons and the right-handed neutrinos Uipt, while ordinary 
multiplication between vectors and matrices is implied Q. 

To obtain a phenomenologically relevant model in this minimal setup, at least 3 
singlet heavy Majorana neutrinos 1^1,2,3 i? are needed and these have to be nearly degenerate 
in mass [13]. To ensure the latter, we assume that to leading order, the singlet Majorana 
sector is SO (3) symmetric, i.e. 

mju = mjvlg + AniA/ , (2.2) 

where I3 is the 3x3 identity matrix and Aitia,/ is a general S0(3)-breaking matrix induced 
by RG effects. As we will discuss below, compatibility with the observed light neutrino 
masses and mixings requires that {AmM)ij/fnN ~ 10"'', for electroweak-mass Majorana 
neutrinos, i.e. for itln ^ 0.1-1 TeV. We will explicitly demonstrate, how such an S0(3)- 
breaking matrix Ahim, of the required order, can be generated radiatively via the RG 
evolution of the right-handed neutrino mass matrix niM from the GUT scale Mx ~ 2 x 
10^^ GeV to the mass scale of the right-handed neutrinos ttin. 

To one-loop order, the RG equations governing the 3x3 matrices of the neutrino 
Yukawa couplings h'^, the charged-lepton Yukawa couplings and the singlet Majorana 

* Occasionally, we will also denote the individual lepton numbers with Le p.,-, but hopefully the precise 
meaning of ^e,/j,r can be easily inferred from the context, without causing confusion. 
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Le,eR Lf^,flR Lr,TR 


1^1 


1^2 ± il^3 




+ 1 +1 





±1 




0+1 









Table 1: Flavour charge assignments for the breaking S0(3) — )■ ^{l)^ x U(l)^ . 



neutrino masses m^/ are given by [25 



dt 

dY/ 

It 

dmM 
dt 



167r2 
1 

167r2 



-g2 + ^9,-T]h 



167r2 



{h^^h") uiM + mM(h"^h'^*) 



(2.3) 
(2.4) 
(2.5) 



where t = ln(Mx//i) and fi is the RG evolution scale. Moreover, gi and g2 are the gauge 
couplings of the U(l)y and SU(2)i gauge groups, respectively, and T is the shorthand 
notation for the trace 



T = Tr(^3h"h"t + 3h'^h^t + h'^h"^ + h^h^t) . 



(2.6) 



Here, h" and h*^ are the 3x3 matrices for the up- and down-type Yukawa couplings, respec- 
tively. At the GUT scale Mx, we impose the universal boundary condition: mM{Mx) = 
m^vls. The corresponding boundary values for the neutrino Yukawa couplings depend 
on the particular model of RiL, which we now discuss in detail. 



2.1 Resonant r-Leptogenesis 

In the physical charged-lepton mass basis, the SO (3) symmetry imposed on the singlet 
Majorana sector at the GUT scale Mx gets explicitly broken by a set of neutrino Yukawa 
couplings to the subgroup of lepton symmetries: U(l)j^ x U(l)^ . The flavour charge 
assignments that give rise to such a breaking are presented in Table [1] 

As a consequence of the U(l)^^^ xU(l)^ symmetry, the neutrino Yukawa coupling 
matrix takes on the general form: 



/ ae"^^/^ 

6e-*^/^ 6e*^/^ 







(2.7) 



/ 
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and Sh'^ vanishes, if the symmetry is exact. In this symmetric hmit, the hght neutrinos 
remain massless to all orders in perturbation theory, whilst a and b are free unconstrained 
complex parameters. The phases accompanying these parameters in (12.71) are simply chosen 
for convenience to maximize the lepton asymmetry in leptogenesis (see Section 14. ip when 
a and b are real. In order to give masses to the light neutrinos, the following minimal 
departure from the flavour symmetric limit is considered: 



( 






-i(7r/4-7i) 






31(71/4-72) 



where |e 



e,/i,T I ) 



1^1,2 ^ Wl, \b\ and the phases 71^2 are unrestricted. A more precise determi- 
nation of the range of parameter values can be obtained from the low-energy neutrino data 
and successful leptogenesis. 

It is important to notice that the flavour structure of the neutrino Yukawa cou- 
plings h" is preserved through the RG evolution, as long as ^h'^ remains a small perturba- 
tion. In detail, RG effects violate the S0(3)-invariant form of the Majorana mass matrix 
mM{Mx) = rriiyl^ by the 3-by-3 matrix Auim- In the leading-log approximation, the 
S0(3)-breaking matrix Auim reads: 



AniM 
/ 

1 

V2 



"87r2 



In 



M 



X 



TUN 



Re [h^^MxWiMx 



rriN 



In 



M 



X 



TUN 



(2.9) 



+ + c 



be, 



erKi sm 7i 



-^(aee + be^) + e^^i sin 71 

0^ + 62 + 



-^{aee + bt^ 

K1K2 sin (71 + 72 
2 



e^K2 sin 72 ^ 



y:^(aee + &e^) + e^K2 sin72 K1K2 sin(7i + 72) a^ + b'^ + nl 



J 



h^(mjv) 



I3 + 



I3 + 



In 

167r2 



ln(Mx 

mjv 



with 71^2 = 7i,2 + f • Correspondingly, the neutrino and charged-lepton Yukawa couplings 
modify via RG running from Mx to mjv as follows: 

fU-Sa"^ -Sab 

—Sab 
\ 
/U' + 3a2 

Sab 
\ 



167r2 



U-Sb^ 


Sab 

u' + 352 





\ 








h-'iMx), (2.10) 



h\Mx) , (2.11) 



U'-lhlJ 



where U and U' stand for the shorthand expressions 

U = -q'i + -q't - Shi - 3/i? - /i^ 

u' = 



9 o So 
^92 + 



2a^ 



9 



15 



1^2 + -j9l 



t 

Shi - Shi - hi - 2a^ 



2b' 



2b' 



(2.12) 
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Observe that the universal contributions U and U' are dominated by the top-quark Yukawa 
couphng ht and are approximately equal, i.e. ?7 ~ f/' ~ — 3. For the models of interest to 
us, we have a, b, h^. ~ 10~^ <^ hf, so the RG effects give rise to an overall rescaling of the 
charged and neutrino Yukawa couplings, and h'^. Without loss of generality, we may 
assume that the charged lepton Yukawa-coupling matrix is positive and diagonal at the 
scale rriN, i-e. close to the electroweak scale. Given the form invariance of under RG 
effects, we may therefore define all input parameters at the right-handed neutrino mass 
scale rriTv, where the matching with the light neutrino data is performed. 

We may now determine the Yukawa parameters (a, b, eg, e^, e,-), in terms of the light- 
neutrino mass matrix m'^ in the positive and diagonal charged lepton Yukawa basis. To do 
so, we first notice that the chosen symmetry U(l)^^^^^ x U(l)j^^ is sufficient to ensure the 
vanishing of the light neutrino mass matrix m'^. In fact, if it is an exact symmetry of the 
theory, the light neutrino mass matrix will vanish to all orders in perturbation theory [26] . 
To leading order in the symmetry-breaking parameters ArriM and ^h'^, the tree- level light 
neutrino mass matrix m'^ is given by 



m 



( 



2mjv 



h'^AniMh 



2m N 



^ab - e,e 



-ab 



Amjv ^2 _ g2 

JTtjv M 



(2.13) 



-4/ 



where v = 2Mw / Qw = 245 GeV is the vacuum expectation value of the SM Iliggs field $. 
5 the las' 

where Atun stands for the expression 



In deriving the last equation in fl2.13p . we have also assumed that \/ ^^^ ^1,2 ^ ^e,/^,r, 



AiriN = 2(AmM)23 + i 



33 



fAniM' 



22 



friN 
87r2 



In 



M 



X 



niN 



2kiK2 sin(7i + 72) + i{K\ ^ '^i) 



(2.14) 



As a consequence of the flavour symmetry U(l)j^^_^^^ x U(l)^^, the symmetry-violating 
parameters ee,^,T and K1.2 enter the tree-level light neutrino mass matrix in f l2.13p 
quadratically. This in turn implies that for electroweak-scale heavy neutrinos ^ v, 
the symmetry-breaking Yukawa couplings 6h^ in (12. 8 p need not be much smaller than the 
electron Yukawa coupling he ~ 10~^. Moreover, one should observe that only a particular 
combination of SO(3)-violating terms ( AmM)ij appears in through ArriN. Nevertheless, 
for electroweak-scale heavy neutrinos with mass differences |AmAr|/mjv ~ 10~^, one should 
have \a\, \b\ < 10^^ to avoid getting too large light neutrino masses much above 0.5 eV. 
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Given the analytic form f l2.13p of the hght neutrino mass matrix, we may directly 
compute the neutrino Yukawa couplings (a, b, eg, e^, e,-), as functions of niN, the phenomeno- 
logically constrained neutrino mass matrix m'^ and the symmetry-breaking parameters k,i^2 
and 71,2: 



2m 



N 



2m 



N 



V 

2m 



\n{Mx/mN) 
ln(Mx/mjv) 

\2 



im' 



13) 



m: 



33 



(m: 



23) 



33 



N 



m' 



33 



m 



11 



m 



22 



(m 



13) 



m' 



33 



m: 



23) 



m' 



33 



2kiK2 sin(7i + 72) + i{K 
2kiK2 sin(7i + 72) + i{K 



(2.15) 



Since the approximate light-neutrino mass matrix in (I2.13P has rank 2, the lightest 
neutrino mass eigenstate vi will be massless in this approximation. The relations given 
in (]2.15p will be used to obtain numerical estimates of the BAU, in terms of m^ and 
the symmetry-breaking parameters k,i^2 and 71^2, for both normal and inverted hierarchy 
scenarios of light neutrinos. In the following, we discuss the two remaining flavour variants 
of RL: R/iL and ReL. 



2.2 Resonant /x-Leptogenesis 



The flavour scenario of R/iL gets realized, once the GUT-scale SO (3) symmetry gets bro- 
ken to U(l)^ X U(l)^ . The flavour charge assignments related to this breaking are 
presented in Table |2j As a consequence, the neutrino Yukawa coupling takes on the 
form: 

t/4 \ 



/ ae-^^/^ ae*^/ 




where the subdominant neutrino Yukawa-coupling matrix 



(2.16) 











\ er 







(2.17) 



breaks minimally the U(l)^^_,_^^ x U(l)^ flavour symmetry. 
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Le,eR Lfj_,flFl Lr,Tf> 


1^1 


1^2 ± il^3 




+1 +1 





±1 




0+10 









Table 2: Flavour charge assignments for the breaking 5*0(3) — t- f/(l)Le+LT x f^(l)L^- 



As in the RrL case, we assume that the Major ana-mass matrix mj\,/ is proportional 
to I3 at the GUT scale Mx and gets radiatively broken via RG effects at the heavy Majorana 
neutrino scale rriN. Taking into account both symmetry-breaking terms Am^ and 6h^, 
the light-neutrino mass matrix acquires an analogous form in R/xL: 



m 



2m 



N 



_ ^ 



■ab 



— e 



— e 



■ab — 



mjv 

AniK T 2 



(2.18) 



where Am^ is given by (I2.14p . Correspondingly, the neutrino Yukawa coupling parameters 
(a, 6, ee, e^, e,-) may be analogously expressed, in terms of m'^, mAr, ^1^2 and 71,2 as follows: 



,2 



2m 



TV 



t;2 ln(Mx/m7v) 



2m 



N 



t;2 ln(Mx/m7v) 



2mj^ {m 



V \2 

12) 



m: 



22 



2mAf ^ 
,2 "^22 ' 



2m AT (m23)^ 



m 



11 



m 



33 



m: 



m: 



22 
•\2 

23/ 



m: 



22 



2kiK2 sin(7i + 72) + i{K 
2kiK2 sin(7i + 72) + i{K 



1 -1 



(2.19) 



m; 



■22 



2.3 Resonant e-Leptogenesis 



A third possible flavour scenario pertinent to ReL is realized by the symmetry-breaking 
pattern S0(3) — )■ \J{l)j^^_^_j^^ x U(l)^_^, where the flavour charge assignments are given in 
Table 131 In ReL, the neutrino Yukawa coupling has the structure: 



/O 







\ 



ae''''/^ ae*^/^ 
\ 6e-*^/^ be^""/^ J 



+ 



(2.20) 
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Le,eR Lfj_,flFl Lr,Tf> 


1^1 


1^2 ± il^3 




+1 +1 





±1 




+10 









Table 3: Flavour charge assignments for the breaking S0{3) — t- U{1)l^j^-l^ x U{1)l^. 
and the breaking terms Sh^ of the U(l)^ x U(l)^^ flavour symmetry are 













(2.21) 



In close analogy with the previous two scenarios of RrL and R/xL, the light neutrino 
mass matrix in ReL is given by 



m 



2m N 



\ 



(2.22) 



where AniN retains its analytic form of fl2.14p . By analogy, we may derive in ReL the rela- 
tions of the neutrino Yukawa coupling parameters (a, b, ee, e^, e,-) to m'^, m^, ki^2 and 71,2- 
These are given by 



2m 



N 



t;2 \n{Mx/mN) 
2m M Stt^ 



(m 



- 



V \2 
12) 



m' 



ln(Mx/mjv) 
2mN ^ 

2mN {m'i^f 



m^ - ^"^^3)' 

''''33 ~ 



m 



2kiK2 sin(7i + 72) + i{Kl - kI) 



2kiK2 sin(7i + 72) + i{Kl - kI) 



(2.23) 



9 u ' 

2m N ("^13)^ 



3 Low Energy Observables 

Low-energy neutrino data provide indisputable evidence for neutrino oscillations. We use 
these data to determine and constrain the fundamental parameters of the theory. We also 
present the full range of predictions for the half hfe of neutrinoless double beta decay of a 
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heavy nucleus, within the different ffavour scenarios of RL. In the same context, we present 
analytic results and estimates for LFV observables, such as /x — )■ 67 and /x — )■ e conversion 
in nuclei. 



3.1 Light Neutrino Oscillation Data 

The interpretation of the experimental results on solar neutrinos suggests that z/g — )■ z/^^,- 
oscillations are mainly driven by the mass squared difference Am^2 = "^^2 ~ "^^i' while 
the corresponding experimental data on atmospheric neutrinos are interpreted by z/^ — z/^ 
oscillations driven by Am^^ = m^^— m^^, within a minimal scheme of three active neutrinos. 
For the present analysis, we use the global fits performed in [27] , even though global fits of 
other groups give compatible results [28]. The best fit values for the light neutrino masses 
and mixings, including their uncertainties at the 2a level, are given by 



sinHu = 0.32+°;°!, sin^^ss = 0.50^°;^^, sin^ ^13 = 0.0071°:°^?, 
Am?2 = (7.6t°:^) X 10-^ eV^ , \Aml,\ = (2.4+°:3) x IQ-^ eV^ (3.1) 



where the sign of Am^g = Am^2 + Am^g remains still undetermined, corresponding to the 
cases of the so-called normal (Am^g > 0) and inverted (Am^g < 0) neutrino mass hierarchy, 
respectively. As outlined in Section [21 we use the oscillation parameters as input to reduce 
the number of free parameters in the neutrino Yukawa coupling matrix h'^. 



Models that include Majorana neutrinos violate the L-number and so can give rise to neutri- 
noless double beta decay (Oz^/3/3) of a heavy nucleus where two single /3 decays [29H3T] 
can occur simultaneously in one nucleus, + 2e~ . The measurement of the half 

life of this decay provides further information on the structure of the light neutrino mass 
matrix m'^. The half life Ty^^ for a Oz//3/9 decay mediated by light Majorana neutrinos is 
given by 



where (m) denotes the effective Majorana neutrino mass, me is the electron mass and 
M oi^fsiB and Gqi are the nuclear matrix element and the phase space factor of the decay, 
respectively. The effective neutrino mass (m) is given by the entry {11} = {ee} of the 



3.2 Decay 




(3.2) 
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light neutrino mass matrix m*^, which can be expressed as 



(m) = 



u\2 



1=1 



(3.3) 




where If^ is the PMNS neutrino-mixing matrix [32j. As described in Section O R£L models 
realise a light neutrino mass spectrum with a vanishing lightest neutrino mass, mj,^ = 0, 
and either a normal or inverted mass hierarchy. The prediction for the effective Majorana 
neutrino mass for these two scenarios is given by 

. , _+^-^) X 10-3 Am?, > 
m) = r ■ (3.4) 

+=^-°) X 10-2 Amj, < 

The above prediction takes into account the uncertainty of the observed oscillation param- 
eters at the 2a level and the variation of the unknown Dirac and Majorana phases. These 
predictions are to be compared to the experimental bound T^!^^ > 1.9 x 10^^ years in 
the isotope ™Ge reported by the Heidelberg-Moscow collaboration [33], implying an upper 
limit on (m) in the range: 

(m)exp < (0.3- 0.6) eV. (3.5) 

Here, the main uncertainty is due to the choice for the nuclear matrix element that occurs 
in (13. 2p . Future Oi//3/3-decay experiments are expected to probe (m) to sensitivities of order 
10-2 [19] and so fall within the range given in (13. 4 p to validate the mass scenario of inverted 
hierarchy. 



3.3 h Z27 

Heavy Majorana- neutrino loop effects may induce sizeable LFV couplings to the photon and 
the Z boson. These couplings give rise to LFV decays, such as yU — )■ 67 [Mj, /i — )■ eee [35] 
and /i — )■ e conversion in nuclei. The strength of LFV is controlled by the effective coupling 
matrix 

f^W. = ^ (h'^h'^^)y. , (3.6) 

which governs the flavour transition between the charged leptons li^2 = e, /x, r in LFV 
processes. The LFV decay li — > Z27, ^1 G {/i, t}, with I2 G {/x,e}, I2 7^ h, whose branching 
fraction is given by 
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In the above, Ti^ is the decay width of lepton li and GH^^ is a composite form-factor given 
by [35] 

G\^'^ = -Qi^i^G, (^) , (3.8) 

with 

= - 4(1 -x)3 -2(1^' ■ ^ 
Given that the experimentally measured muon and tau decay widths are = 2.997 x 
10~^^ GeV and = 2.158 x 10~^^ GeV [36], the LEV branching ratios can be expressed 
as 

Biii^ej) ^ 8.0 ■ 10-^ X (7 f , (3.10) 

B{r^hi) ^ 1.5 ■ 10-^ X f \^ri,W l2 = e,fi. (3.11) 

Here, we defined g{x) = 4G^(x^), which possesses the limits (7 —t- 1 for ^ and 
(7 — )■ 1/16 for m^v = ""^VK- Our theoretical predictions will be contrasted to the current 
experimental upper limits [36] 
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5exp(/i ^67) < 1.2 X 10" 

Sexp(r ^ /i7) < 6.8 X 10-« , (3.12) 
5exp(r^e7) < l.lxlO-^ 

and the expected sensitivity of the MEG experiment [20] , 

SMEG(/i^e7) ~ 10-^^ (3.13) 

In RiL models, only two of the right-handed neutrinos, z/2_r and u^r, have appreciable e- and 
/X- Yukawa couplings, a,b ^ 10^^, ki^2 = 10^^ — 10"^, and will be relevant to LEV effects. 
Eor example, in the RrL model, the parameters lil/^^jp are, to a good approximation, given 
by 

44 4 

I/-V |2 ^ 27 2 \r-\ |2 ^ / 2 2\72 i /-» |2 / 2 2\ 2 

(3.14) 

Because of the relations (12.151) . a and h are approximately inversely proportional to ki^2, 
a, 6 ~ 3 X IQ^"^ (in the case of a normal light neutrino mass hierarchy), and for a typical 
value of Ki 2 ~ 10~^, a and h are of the order a, 6 ~ 3 x 10^'^. Thus, in the RrL scenario. 
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the following typical values for the branching ratios of the photonic LFV jj and r decays 
are predicted: 

B{fi^e^)KrL ~ 10-^', 5(r^/27)RrL ~ 10-^^ (3.15) 

Consequently, only — > 67) is expected to be within reach of future experiments, 
whereas the LFV r decays are far beyond the realm of detection. By analogy, the pre- 
dictions in the Re(/i)L model are found to be 

5(/i^e7)ReML ~ 10-^', 5(r^/27)ReML ~ 10-^'. (3.16) 

Again, i?(r — )• fx^y) and B{t ^ 67) are expected to be below current and future experi- 
mental sensitivities for parameter choices compatible with successful leptogenesis. While 
the prediction for i?(/i — > 67) is far below the expected MEG sensitivity, the LFV /i — ?• e 
transition rate might still be high enough to be testable at future experiments measuring 
/i — )■ e conversion in nuclei (see next section). Because of the inverse proportionality be- 
tween a, h and ^1,2 the prediction for i?(/i — 67) is essentially independent of the choice 
for K12. In Section [5l we present detailed numerical results for LFV /i- and r-decays that 
confirm the validity of these simple estimates. 



3.4 Coherent — )- e Conversion in Nuclei 



One of the most sensitive probes of LFV is the coherent conversion of /i — )■ e in nuclei 
The /i — )■ e conversion rate in a nucleus with nucleon numbers (A^, Z) is given by 

V[ii{N,Z)^e{N,Z)] 



B,eiN,Z) = 



3 4 5 



74 

\F(-m 



2m2 



\Q 



w 



(3.17) 



r[/i (A^, Z) capture] IGTr^m^Fcapt Z 

where = 1/137 is the electromagnetic fine structure constant, Z^g is the effective atomic 
number and Fcapt is the muon nuclear capture rate. For 22^1, experimental measurements 
give Zeff ^ 17.6 [40j and F^apt ~ 1-705 x 10"^^ GeV jS]. Moreover, \F{-ml)\ ^ 0.54 is 
the nuclear form factor Il2l and 



Q 



w 



Vu{2Z + N) + Vd{Z + 2N) 



(3.18) 



is the weak matrix element, where 



Vu 



1+3^ + 
'1 _ 1 2 

4 a'"'^ 



8 9 

3 _ 4 2 
8^9''"' 



In!^ 



m 



w 



In 



m 



N 



m 



w. 



1 fl 



Me + 2 



4 + 3< 



1 fl 



4 3''"' 



m 



N 



m 



m 



w 



N 



m 



/ fie 



W 



(3.19) 
(3.20) 
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In the case of 22^1, -B^e(26, 22) is then approximately related to B{fi e-y) through 

S^e(26,22) ^ 10-1 X B{fi ^ 67), (3.21) 

for a right-handed neutrino mass scale itin ~ 100 GeV. On the experimental side, the 
strongest upper bound is obtained from data on /i — )■ e conversion in 22"^^ PS]) 

5;^^P(26,22) < 4.3 X 10-1^ (3.22) 

Using the relation (13.211) . we observe that this sensitivity is comparable to the current bound 
on i?(yU — 7- 67). However, the proposed COMET and mu2e experiments, measuring — )■ e 
conversion in 13AI, are expected to be sensitive to conversion rates of order 10"^^ [441145] , 
thereby improving the sensitivity compared to the current limit by four orders of magnitude. 



4 Leptogenesis 

In this section, we briefly review the central results of the field-theoretic formalism for RL 
developed in [7l[8] which will be used in our analysis. We then set up the BEs and present 
approximate solutions for the kinematic regime of interest to us. Finally, we clarify the 
necessity of introducing at least three right-handed neutrinos into the theory, in order to 
explain the BAU and obtain testable rates of LEV, such as — i- 67) ~ lO^i^-lO^^^. 

4.1 Leptonic Asymmetries 

Within the framework of leptogenesis, a net non-zero leptonic asymmetry results from 
the CP-violating decays of the heavy Majorana neutrinos Na into the left-handed charged 
leptons l2 and light neutrinos uil. Consequently, we have to calculate the partial decay 
width of the heavy Majorana neutrino into a particular lepton flavour I, 

r^i = r{N^ + W^) + T{N^ ^ + Z,H). (4.1) 

Eor temperatures above the electroweak phase transition, the SM VEV vanishes and only 
the would-be Goldstone and Higgs modes will predominantly contribute to Tai- 

In RL models, resumming the absorptive parts of the heavy Majorana- neutrino self- 
energy transitions N13 — )■ Na plays an important role in the computation of Tai [HIH]. In 
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order to take this resummation consistently into account, we first introduce the lepton- 
fiavour dependent absorptive transition amphtude [36] 



Al,{hn ^ , (4.2) 

which represents the contribution of a single charged lepton and light neutrino flavour / 
running in the loop. Summing over all flavours /, we then get the total transition amplitude 

A^^ihn ^ J2 ^^U^^) = iJ"' ■ (4-3) 

l=e,fi,T 

Note that the diagonal transition amplitude Aaa is related to the tree-level decay width of 
the heavy Majorana neutrino Na through: = 2mN^Aaa(h'^)- 

Since all charged-lepton Yukawa couplings will be in thermal equilibrium in the low- 
scale leptogenesis scenarios of our interest, we consider the weak basis in which the matrices 
and uim are both diagonal and positive. To account for unstable-particle mixing effects 
between the 3 heavy Majorana neutrinos, we follow (TJIH] and define the resummed effective 
Yukawa couphngs h^^ and their CP-conjugate ones h^^ related to the vertices L^Na and 
L'-^^*Na, respectively. The resummed neutrino Yukawa couplings h^^ are given by [8|H6] 



iha^ = K - ^ E M (4.4) 

/3,7=1 

^ rriaimaAa^ + mpAp.c) - iRa-y [maA^i3{maAa-y + m^A^a) + rni3Aii^{maA^a + '"^7^07)] 
m%-nip + 2imlAf^i3 + 2ilmRa^ (m^ + m^m-^ReA^^) 

where \eai3'y\ is the modulus of the usual Levi-Civita anti-symmetric tensor (£123 = 1), 
= "^L' ^<^P = ^afiih") and 

Ral3 = — n 2 i"o- — • ('^•^) 

The respective CP-conjugate effective Yukawa couplings h^^ {h'^) are obtained from (14. 4p 
by replacing the tree- level couplings hj'^ with their complex conjugates hj'^: 

(h^) = ihaii^n ■ (4.6) 

In our calculations, we neglect the 1-loop corrections to the proper vertices L^Na, whose 
absorptive parts are numerically insignificant in RL. In terms of the absorptive transition 
amplitudes A'-^^ given in (14. 2p . the partial decay widths P^^ and their CP-conjugates P^^ 
may now be expressed in the compact form: 

^ai = ^iv.AL(r) , P^, = m^„AL(r^) , (4.7) 
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where we have exphcitly indicated the dependence of the absorptive transition amphtudes 
on and h'^'" . Given the analytic expressions (14. 7p . it is straightforward to compute the 
leptonic asymmetries for each individual lepton flavour: 

S^, ^ tjj = 1^^"' 1^^°' ^ . (4.8) 

l=e,^J.,T ^ ' 

The analytic results for the leptonic asymmetries 8cd simplify considerably in the 2-heavy 
neutrino mixing limit, in which defined in (14. 5p is set to zero. In this limit, 5q/ are 
given by [ZIIH] 

^ Im[ (h-|hr,) (h-th-).,] {ml^ - m^^) m^, rj;) 

where a, /3 = 1,2 and rj^^ is the tree-level decay width of the heavy Majorana neutrino Na-, 
given after (14. 3p . Based on the simplified expression (14. 9p . the following two conditions for 
having resonantly enhanced leptonic asymmetries bod. ~ C'(l) may be derived [7]: 

(i) |m;v, -m^J ^ , (4.10) 



|lm[(h-|hr,)(h-th-).,] 

(h^^h^),, (h'^th-^ 



- 1 . (4.11) 

Note that the first resonant condition (i) is exactly met, when the unitarity limit on the re- 
summed heavy-neutrino propagator gets saturated |4^, i.e. when the regulating expression 
in (USD, 



12 2 1 -nW 



f = ' "° '2^ " < 1 (A 12) 



takes its maximal possible value: /reg = 1. Within our RL scenarios, the first condition 
in (14.1 op is naturally fulfilled as the heavy-neutrino mass splittings are generated via RG 
effects and are of the required order. The second condition is crucial as well and controls 
the size of the leptonic asymmetries. As we will see below, the condition (ii) in (I4.10p has 
a non-trivial impact on approximate L-conserving RL models. 



4.2 Comparison with Other Methods 

It is now worth commenting on some of the attempts made in the literature P81450] to 
calculate the resonant part of the leptonic asymmetries dai- Their results differ by the way 
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in which the singularity rriN^ — )■ rriN^ occurring in the denominator of the second fraction 
in f l4.9p gets regulated, when heavy-neutrino width effects are taken into consideration. 
Specifically, the various approaches differ in their derivations for the analytic form of /reg 
given in (14.121) . For instance, the authors of [H] use a perturbative quantum- mechanical 
approach to obtain a regulator of the form: 

_ Am^rg;^/2 

wheremTv = ^{171^^+171^^) and AmAr = |mAr^— mTVal- It is easy to observe that for scenarios, 
for which Re(h'^"''h'^)i2 = 0, but Re(hJ'/hJ^2 ) 7^ 0; the unitarity upper bound given in (14.121) 
gets violated, in the degenerate heavy- neutrino mass limit rriN^ ttln^. In particular, in the 
same limit, the individual lepton-flavour asymmetries 5ii (with I = e, fi, r) become singular. 
Although this singularity disappears when the lepton-flavour sum 6ni^2 = X]/=e t ^i,2Z is 
taken, the regulator (I4.13P will still be inapplicable to lepton-flavour RL scenarios, for which 
AmN/rriN ~ Re(h'^th'^)^2/(167r) < TN.^JmN- 

Based on a modified version of the field-theoretic approach introduced in [71|8], the 
authors of |l9l|50] obtain a different regulating expression for the leptonic asymmetry ^i;: 

12 21 -p(o) 

■I^^S / n n ^2 / „(0) p(0)n2 • l^i'-L^J 



It is not difficult to observe that /reg diverges as {rriN-^ — rriN^)'^ in RL scenarios, for which 
(h''^h.'^)u = (h'^'''h'')22, even though one could still have (h^/hj'^ ) 7^ (h.'j^2^'^2) G^ch single 
lepton flavour /. For instance, such a situation can naturally occur in approximate lepton- 
number conserving models of RL discussed recently in [SI1E2]. Evidently, the prediction 
for the leptonic asymmetries ^tvi 2 such scenarios may get overestimated by many orders 
of magnitude. 

In order to illustrate this last point, let us consider a one-generation inverse seesaw 
model [53] with two singlet neutrinos of opposite lepton number, and Sl- The neutrino 
sector of this model is described by the lepton-number-conserving Lagrangian 

- inverse = ^ {^L, U^) ^ ^ ^ ^ + ^i? (^L, h) ^l^R + H.C. (4.15) 

Without loss of generality, the kinematic parameters M and hn can be rephased to become 
real. Following closely the discussion in [7], we introduce into the Lagrangian (I4.15P the 
lepton-number violating operators | ^rUj^Uj^, ^ fj,L S^S^ and hi {i^l, Il)^S1. In order to 
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minimally break both the lepton number and CP, it was shown in [7] that at least two of the 
aforementioned AL = 2 operators are needed. In fact, this result is a direct consequence 
of the Nanopoulos- Weinberg (NW) no-go theorem [5^. The NW theorem states that no 
net baryon asymmetry can be generated by a single B- and CP-violating operator to all 
orders in perturbation theory. 

It is interesting to provide an estimate of the lepton asymmetry obtained within a 
simple model of approximate lepton-number conservation, where fii = —ilJ-R = A* and 
hi = 0. For these parameters, we adopt the same ballpark of values as in [51]: ttln ^ 
M = 1 TeV, hn = 3 X 10"^ and Atun ^ /i/v^ = 2 x 10"^° M. With these input 
parameters, we may estimate that P^^,^ 2 x lO^^mTv, leading to AmTv/Pj^j^ ~ 10~^. 
The first fraction containing the CP-violating phases in fl4.9p is rather suppressed, of order 
/x/M ~ AmTv/'^AT ~ 10"-*^", whilst /reg ~ 2Am7v/PAri 2 ~ 2 X 10~^ within our resummation 
approach. This gives rise to dismally small lepton asymmetries: 

Sii ^ 621 ^ 10-'' . (4.16) 

Notice that the lepton asymmetries 61^21 are proportional to /i^ in agreement with the NW 
theorem and, as they should, both vanish identically in the L-conserving limit — 0. 
As we will see in the next section, lepton asymmetries of order 10^^^ will fall short by at 
least 7 orders of magnitude to explain the BAU. Instead, had one used the regulator /reg 
in fl4.14p . one would have obtained the enormous enhancement: f^eg ~ r7Vi.2//^ ~ 10^ ^ 1, 
leading to lepton asymmetries of order [51\: 61^21 ~ ^Ni,2/™n ~ 10~^. However, this 
result is independent of the L-conserving parameter = fi^ = —ifiR, and taking the L- 
conserving limit /i — )■ 0, one obtains a non-zero lepton asymmetry, which clearly signifies 
an erroneous result. The above exercise shows that one has to go beyond the two-heavy 
neutrino mixing framework and consider non-trivial lepton flavour effects [13], in order to 
obtain a phenomenologically relevant model that predicts testable rates for LFV. 

4.3 Baryon Asymmetry 

In this section we present the relevant Boltzmann equations (BEs) which will be used 
to evaluate the heavy-Majorana- neutrino-, the lepton- and the baryon-number densities, 

^Afi,2,3 = ''^A^i,2,3/''^7) ''^ie.M.r = ""-Le.M.T /''^T T] B = n^/^^, uormahscd to the photon number 
density n^. In our computations, we include the dominant collision terms related to the 
1 — )■ 2 decays of the right-handed neutrinos and to resonant 2 — )■ 2 lepton scatterings that 
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describe AL = 2 and AL = transitions. We neglect chemical potential contributions 
from the right-handed charged leptons and quarks. 

In order to appropriately take into account the flavour effects on the BEs for r/Ar^ 2,3 
and ?7l^^^, we closely follow the approach and the notation established in More 
explicitly, the BEs are found to be 



dr]Li _ £ 

dz n^H 



eq 

'/ 

3 



7^1 , (4.17) 



3 "^L, (^7ip$t - 7l'<i. j 



k=e,iJ,,T 



(4.18) 



where a = 1, 2, 3 and I = e, fi, t. In addition, ~ 17xm^/Mp is the Hubble parameter 
at T = rriN, where Mp = 1.2 x 10^^ GeV is the Planck mass. The T-dependence of the 
BEs ( I4.17P and (I4.18p is expressed by virtue of the dimensionless parameter 

. = ^ . (4.19) 
in terms of which the photon number density is given by 

. ^ C(3) ^ ^ . (4.20) 

TT^ TT^ Z-^ 

where C(3) ~ 1.202 is Apery's constant. Finally, 77^ in fl4.17p and fl4.18p is the equilibrium 
number density of the heavy Majorana neutrino N^, normalised to the number density of 
photons, i.e. 

^ \z'K, (z) , (4.21) 
where Kn{x) is the nth-order modified Bessel function j55j . 

The BEs (14.1 7p and (I4.18P include the collision terms for the decays — > Li^, as 
well as the AL = 0,2 resonant scatterings: Lfc$ o L/$ and Lfe$ o Lf^\ which are 
defined as [H] 



j^^fi ^ 7(^.$ ^ m + 7(^^$^ ^ Lf<^^) , (4.22) 
7Xt = ^iLk<!> ^ LP$t) + 7(L^$t ^ L,$) . 
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Since we are only interested in the resonant part of the above 2 — 2 scatterings, we make use 
of the narrow width approximation (NWA) for the resummed heavy-neutrino propagators. 
Thus, at the amphtude squared level, we employ the NWA for the complex-conjugate 
product of Breit-Wigner propagators [56] : 



1 



Sis - rn%J + 5{s - m%) 



m 



No, 



(4.23) 



Note that the NWA relation fl4.23p is valid for any range of parameter values for Fat-^ j 3 
and long as Fat^ 2 3 ^ ^a^i 2 3- This last condition is naturally fulfilled within the 

R£L models. 

The collision term pertinent to the heavy-Majorana-neutrino decays 7^1 is given by 



7fi 



m 



N 



KAz)T 



(4.24) 



whilst the corresponding collision terms for the AL = 0,2 lepton-flavour transitions are 
calculated by means of the NWA fl4.23p to be 



a,B=l 



^ y^la "-ka "-ip "-kp + "Za "-ka "-113 "-kji 

(r^F).. + (F^^r^).. + (r^F),, + (F^^r^),, 



X 1 - 2z 



Fat, + Fat^ 



-1 



(4.25) 



7LP<I.t 



^ y^la ^ka "-ip "-kp, + "-lex "-ka ^Z/3 



a3=\ 



X 1 - 2i 



h h 

Fa^c + ^Nfi 



I act 
-1 



+ (h' 



Y-vC]rr-vC. 



.rr-vC]Y-vC. 



+ (h-h^),, + (h^-'h^ 

(4.26) 



The scattering collision terms (14.25^ and f l4.26p contain all contributions from the resonant 
exchange of right-handed neutrinos in the NWA, including contributions from the so-called 
real intermediate states (RISs) [57J. The latter are obtained when only the diagonal terms 
a = (3 are taken in the summation over the heavy Majorana states N^^p- 
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Separating the diagonal a = P RIS contributions from the off-diagonal a ^ (3 terms 
in the sum, we may rewrite the BE (14.181) in the form [14j 



dz 



3 



1- a=l 



k=e,fi,T 



k=e,^,T 



1- a=l 



1 Li<l> 



(4.27) 



where 7'y = 1y ~ (7y)ris denote the RIS-subtracted collision terms and Bai are the 



branching ratios 



al 



-V I 2 



+ h 



la 



X] ^ ak + ^'ak) 
k=e,^,T 



(h h)_ + (h ) 



(4.28) 



The collision terms proportional to rj^^ and to SaiSak on the RHS of (14.271) turn out to be nu- 
merically negligible for the R^L scenarios under consideration. Instead, the RIS-subtracted 
collision terms proportional to 77^^ in (I4.27P become significant. Their importance in RL 
models was originally raised in [SlIIl] and confirmed most recently in [51j . 

The next step is to include the effect of the {B + L)-violating sphaleron processes [B] . 
For temperatures T larger than the critical temperature ~ 135 GeV of the electroweak 
phase transition, the conversion of the total lepton-to-photon ratio J2i=e firVLi to the 
baryon-to-photon ratio ri% at Tc is given by the relation: 



28 
51 



l=e,fi,T 



(4.29) 



For T < Tc, the so-generated baryon asymmetry 77^ gets diluted by standard photon 
interactions until the recombination epoch, leading to the BAU 



Vb ~ ^Vb- 



(4.30) 



This theoretical prediction can now be compared with the current observational value for 
the baryon-to-photon asymmetry [2], 



(6.20 ±0.15) X 10 



-10 



(4.31) 



It is instructive to derive approximate numerical solutions to the BEs (14.171) and (I4.27p . 
To this end, we re-express the BEs in terms of the out-of-equilibrium deviation parameters 
^Vn^ = VNc/v'n ~ ^ neglect the suppressed 0{S'^i) colhsion terms. For simplicity, we 
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Figure 1: Numerical solutions to BEs i4-32\) and lji4-35\ ) for 6r]N^ = VNi/Vn ~ ^ ^'^^ Vl^j 
respectively, in a RtL model with = 220 GeV, Ki = 10®, Kf^ = 10^, 6r = 10~® and 
Zc = niN/Tc = 1.6. 



initially ignore the RIS subtracted collision terms. But, as we will see below, their inclusion 
is straightforward. With these approximations and simplifications, the BEs may be recast 
into the compact form: 



d6r] 



No, 



dz 



3 



1 + (l - ^a^) S7]No^ 



dz 



z 



(4.32) 



(4.33) 



with Ka = rAf„/[C(3) Hn]- In the kinematic regime z > z" ^ 2Ka the solution to (14.321) 
may well be approximated by 

SvnA^) = {K^z)-\ (4.34) 
independently of the initial conditions (see Figured]). In this regime, the BE (14.331) becomes 



dz 



(4.35) 
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with 5i = X]a=i "^a' ^^'^ ~ ^a=i ^aBai- We may include the numerically significant 
RIS-subtracted collision terms proportional to rji^ in f l4.27p . by rescaling Ki — )• niKi = Kf^, 
where 



= 







(h h + (h h )al3 




2 


[(h h ')n + (h h ')n 






)/3/3 



' ^ 2(7^^^^)-^/=^ , 4 ~ 1-25 In {25Kf) . (4.37) 



2 



For z values in the range Z2 < z < z^^, the solution riL^ may well be approximated by 



ri..(^) - I 7^ . (4.38) 



For 2; > Z3, the lepton number density rj^ freezes out and approaches the constant value 
?7L; = (35i)/(2ii'f''^2;3)§. The general behaviour oir]^ in the different regimes is displayed in 
Figure [TJ 

In this paper we only consider KiL scenarios, for which the washout is strong enough, 
such that the critical temperature Zc = mi^/Tc where the baryon asymmetry rjB decouples 



^The onset of the freeze-out is defined as the position z\ where the relative slope of ()4.35p drops 
below 1, i.e. when \r]']^JriLi \ = '^/'^{^'hY = 1- The solution to this equation can be analytically 
expressed in terms of the Lambert W function, to which Zg in (|4.37p proves to be an excellent interpolating 
approximation over the wide range of values Kf^ « 2 — 10^°. 



a,/3=l ^ 

x(l-2,:!=^^)"\ (4,36) 

In determining the scaling factor k^, we have assumed that r]^ S> in fl4.18p . which is a 
valid approximation within a given R£L scenario under study. Note that if only the diagonal 
a = P terms representing the RIS contributions are considered in the sum, ki reaches its 
maximum value, i.e. k/ = 1 + 0{5f). We also have checked that in the ///-conserving limit 
of the theory, the parameter vanishes, as it should. 

As is illustrated in Figure [H the solution rj^ to fl4.35p exhibits different behaviour 
in the three kinematic regimes, characterized by the specific values of the parameter z = 
ttin/T: 
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from the lepton asymmetries r/L, is situated in the hnear drop-off or constant regime. 
Specifically, we require that 

> 2K-^I\ > 2{Kf)-'/\ (4.39) 

for all heavy neutrino species A^^^ = and lepton flavours / = e, fi, r. As a con- 

sequence, the baryon asymmetry r]B becomes relatively independent of the initial values 
of riLi^ and r]N^. In this case, taking into account all factors in f l4.29p . f l4.30p and f l4.38p . the 
resulting BAU is estimated to be 

28 1 3 ^ Si 

" ~51 27 2^£^^^^Kfmm{z,,zi) 

~ _ 3 . 10-2 J2 f T ■ (4.40) 



=e,tM,T Kf mm 



m^/Te, 1.251n(25K 



We note that the formula (14.40 1) provides a fairly good estimate of the BAU rjB to less 
than 20%, in the applicable regime of approximations given by fl4.39p for Kf^ > 5 for a 
right-handed neutrino mass scale of the order of the EW scale. Hence, to account for the 
observed BAU r]'^^ given in (I4.3ip . lepton asymmetries 61 ^ 10^^ are required. In the 
next section, we present numerical estimates of the BAU for R£L scenarios, based on the 
simplified BEs (I4.32p and fl4.35p . with Ki replaced by Kf^ as defined by means of fl4.36p . 



5 Numerical Results 

In this section we present numerical estimates of the BAU rjB and the low-energy LEV 
observables for fi ^ e and r (e, fi) transitions, based on the analytic results derived in 
Sections [3] and |H Our aim is to delineate the viable parameter space of the R^L models: 
RrL, R/xL and ReL, by considering both cases of a normal and inverted hierarchy for the 
light neutrino mass spectrum. In all the KiL scenarios under study, the lightest neutrino 
is massless. As described in Section [21 we use the neutrino oscillation data to determine 
the theoretical parameters of the light neutrino mass matrix m'^. Specifically, we invert the 
seesaw formula and solve for the neutrino Yukawa couplings a, b and €e,fj,,T- In addition, 
the electroweak-scale flavour structure of the right-handed neutrino mass matrix ium is 
generated from a flavour-universal heavy-neutrino mass matrix mM{Mx) = '^Afls at the 
GUT scale Mx, after taking into account RG- running effects. 

Eor a given light neutrino mass matrix m'^, the solution obtained for a and b is unique 
up to a common sign factor. This sign degeneracy could be eliminated, only if the sign 
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Figure 2: T/ie baryon asymmetry rjB (blue contours) and the LFV observables B{fi ^ cy) 
and i?^e(22Ti) (red contours) as functions of Ki and k,2 in the RtL model with m^v = 
120 GeV, assuming a normal light neutrino mass spectrum. The remaining parameters are 
chosen to be: 71 = Svr/S, 72 = tt/2, 0i = vr, 02 = 0, Re{a) > (left panel); 71 = n/2, 72 = 
Svr/S, 01 = 0, 02 = 0, Re{a) < (right panel). The neutrino oscillation parameters are set 
at their best fit values, with sin^ ^13 = 0.033 at its 2cr upper limit. The blue shaded regions 
denote the parameter space where the baryon asymmetry is larger than the observational 
value r]f^ = 6.2x10-^^. The green shaded areas labelled as '5m'^ = 0.25 ' and '0.10 ' indicate 
the parameter space where the inversion of the light-neutrino mass matrix is violated at the 
25% and 10% level, respectively. 



of Re (a) were known. There is also a similar freedom for an overall sign in the ee,^,r 
parameters, but this turns out to be irrelevant, since it applies to a complete column in 
the neutrino Yukawa matrix and can be rotated away. In addition to the light neutrino 
masses and mixings, the light neutrino mass matrix may also contain the CP-violating 
Dirac phase 5 and the Majorana phases 0i,2- As our ansatz for the neutrino Yukawa 
matrix \\^ uses maximal CP phases for a and 6, the inclusion of the light neutrino CP 
phases is not expected to increase or change significantly our predictions for the baryon 
asymmetry rjB- We therefore consider only the extreme cases of CP parities, 5, 0i^2 = 0, vr, 
as natural choices which can help us to reduce the dimensionality of the parameter space. 
As well as sin^6'i2, sin^6'23, sin^6'i3, Am^2) l^^^isl' ^^S^ (^"^13)' ^ ^^"^ ^^1,2^ we also 
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Figure 3: The same as in Figure\^ hut assuming an inverted light neutrino mass spectrum. 
The remaining parameters are chosen as follows: 71 = Svr/S, 72 = i^/2, 0i = 0, 02 = 0, 
Re{a) > (left panel); 71 = 7r/2, 72 = Svr/S, 0i = 0, 02 = 0, i?e(a) > (right panel). 

have the free parameters ki 2, 71,2, the sign of Re(a) and m^v in our RiL models. Unless 
otherwise stated, we use the best fit values of [27] for the measured neutrino oscillation 
parameters. 

To start with, we show in Figure |2] numerical estimates of the baryon asymmetry tjb 
and the LFV observables B{fi ^ 67) and _B^e(22Ti), as functions of the Yukawa-coupling 
parameters ki and k,2, in a RrL model with mi^ = 120 GeV and a normal hierarchical 
light neutrino mass spectrum. The neutrino mixing angle sin^ 613 is set to its upper exper- 
imental 2(T limit [27]: sin^ ^13 = 0.033. The remaining theoretical parameters are chosen 
to maximise the overlap between successful generation of the required baryon asymme- 
try and high LFV rates. For definiteness, we set 71 = 37r/8, 72 = l/27r (left panel) and 
7i = l/27r, 72 = Sn/S (right panel). The blue shaded areas denote the parameter space 
where the numerically predicted baryon asymmetry rjB is larger than the observational 
value r]^^*^ = 6.2 X 10"^°. These areas should be regarded as representing regions of vi- 
able parameter space for successful leptogenesis, given the freedom of re-adjusting the CP 
phases 71,2 of the Yukawa couplings ki,2- 

The parameters ki^2 are varied within the range 10~^-10~^; smaller values of k,i^2 
would lead to too large values for a,b > 0.1, whereas larger ki^2 values would violate the 
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Figure 4: The branching ratio B[jj, — > 67) as a function of sin^ 613, in RtL scenarios with 
a normal (blue) and inverted (red) light neutrino mass spectra. The model parameters 
for these two cases are as in Figures (left panel) and (left panel), respectively, with 
Ki = 10^^ and K,2 = 10~^. The horizontal solid (dashed) line denotes the current (expected 
future) limit on B{fi — > cy). The vertical solid and dashed lines denote the 2a upper limit 
and the nominal best fit value of sin^Ois, the latter roughly corresponding to the expected 
sensitivity of future oscillation experiments. 



assumed approximation K1.2 ^ a, b, required to invert the seesaw formula. The degree of 
violation of this approximation is indicated by the green shaded areas in Figure [21 These 
green shaded areas denote the parameter space, in which the quantity 

sm?, = =^ ( ."r^y:^':^/ i (5.1) 



I ; I f I ! I 



rriN- \min(|ee 

is bigger than 0.25 (dark green) and 0.10 (hght green), respectively. The quantity 6m'^ is 
a measure that quantifies the accuracy of our analytic approximation for neglecting the 
contribution of the neutrino Yukawa couplings ki^2 in the light neutrino mass matrix m'^ 
[cf. f l2.13p ]. Hence, the values of 6m'^ = 0.10 and 0.25 mean that the inversion of the 
light neutrino mass matrix is accurate at the 10% and 25% level, respectively, due to the 
assumed absence of the K12 terms. It should be stressed that the corresponding parameter 
space is not ruled out; larger values of ki,2 ^ 5 x 10"^ would require a numerical approach 
to invert the light neutrino mass matrix m'^, beyond our analytic approximation. 
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Figure 5: The baryon asymmetry rjB as a function of the ratio Ki/k,2 o RtL scenario, 
where k,ik,2 = 10^® is fixed. All other parameters are as in Figure IB (left panel). The 
solid line corresponds to the full numerical solution of tjb using the three-heavy-neutrino 
mixing formula given in Iji4-4\ l> the dashed line describes the two-heavy-neutrino mixing 
approximation with R^p = 0, and the dotted line is obtained by neglecting the RIS- subtracted 
collision terms in the BE ^.27 ), or equivalently by setting the parameter defined in ^4 -36 ) 
equal to 1 in the BE (J^ 



The specific choice for the phases 71,2 in Figure [2] approximately maximises the numer- 
ically predicted baryon asymmetry r]B compatible with testable LFV decay rates. Specif- 
ically, the areas around (ki,K2) ~ (10-^,10-^) (left panel) and (10-^,10-^) (right panel), 
where cy) ~ 10^^^ can be achieved, are only present for sufficiently asymmetric 

values 7i 7^ 72, and only for specific choices for the remaining discrete parameters in the 
RrL model. As a consequence, the requirement of both a successful generation of the 
baryon asymmetry rjB and potentially observable LFV rates for /i — e transitions puts 
severe constraints on the model parameter space. We note that the dependence of the 
baryon asymmetry rjB on the right-handed neutrinos mass scale m-Ar is weak in the phys- 
ically interesting region of = 100-500 GeV. Finally, the LFV r decays are extremely 
suppressed with B{r — > /27) ~ 10^^^ {I2 = e,fi), and so remain far beyond the realm of 
detection. 

In Figure |3] we display numerical estimates of 77^ and the LFV observables i?(/i — )■ ej) 
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Figure 6: The baryon asymmetry tjb (blue contours) and the branching ratio B(t — )■ 67) 
(red contours) as functions of Hi and 1^2 the RfiL model with = 120 GeV, where a 
normal light neutrino mass spectrum is assumed. The remaining parameters were chosen as 
follows: 7i = Sn/S, 72 = vr/2, 0i = 0, 02 = 0, Re{a) > (left panel); 71 = 7r/2, 72 = Svr/S, 
01 = 0, 02 = IT, Re{a) < (right panel). The neutrino oscillation parameters are set at 
their best fit values, with sin^ ^^13 = 0.033 at its 2a upper limit. The blue shaded regions 
denote the parameter space where the baryon asymmetry is larger than the observational 
value rjB = 6.2 x 10~^°. The green shaded areas labeled as '6m'^ = 0.25 ' and '0.10 ' indicate 
the parameter space where the inversion of the light-neutrino mass matrix is violated at the 
25% and 10% levels, respectively. 

and B^i,{^Ti), as functions of Ki and ^2, in a RrL model with m^r = 120 GeV and 
an inverted hierarchical light neutrino mass spectrum, characterised by Am^g < 0. As 
before, we set the neutrino mixing angle sin^6'i3 to its upper experimental 2a limit [27] : 
sin^ ^13 = 0.033. As can be seen from Figure[3], the regions of the parameter space that yield 
successful baryon asymmetry rjB are smaller than those found in the normal hierarchical 
case for the light neutrino mass spectrum. If such a scenario gets realized in nature, then 
successful leptogenesis implies rates for B^j,e{^Ti) ^ 10^^^, which can still be within reach 
of the projected PRISM experiment [2T] . 

The higher rates for the LEV /i — )■ e transitions in a RrL model with an inverted 
hierarchical light-neutrino mass spectrum may be attributed to the fact that the squared 
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Figure 7: The same as in Figure 0, but for an inverted light neutrino mass spectrum and 
the following choice of parameters: 71 = 3n/8, 72 = tt/2, 0i = 0, 02 = 0, Re{a) < (left 
panel); 71 = n/2, 72 = Sn/S, 0i = 0, 02 = 0, Re{a) > (right panel). 

parameter is proportional to the not yet well determined neutrino-mixing angle sin^ ^13. 
In detail, the branching ratio — )■ 67) is enhanced in the case of an inverted hier- 
archical light neutrino mass spectrum {Amf^ < 0). Instead, in the normal hierarchical 
light-neutrino scenario with Amf^ > 0, — )■ 67) is essentially independent of sin ^13. 
This is demonstrated in Figure HI showing B{fi ^ 67) as a function of sin^6'i3, where 
Ki = 10~^ and K2 = 10~^. In the same figure, the vertical solid and dashed lines denote 
the current best fit value and the expected sensitivity of future experiments for measuring 
sin^6'i3, respectively. The branching ratio B{fi — 67) depends linearly on sin^ 6'i3 in the 
inverted hierarchical light-neutrino scenario where Am^g < 0. Unlike in the normal light- 
neutrino scenario, the predicted value for r^^ in the RrL model with an inverted hierarchical 
light neutrino mass spectrum falls short of explaining the BAU by two orders of magnitude, 
for a potentially observable branching ratio of 67) ~ 10^^'^. 

In Figure [5] we illustrate the importance of including the full three-heavy-neutrino 
mixing in the calculation of the baryon asymmetry. We present a comparison between 
the full calculation based on using the effective Yukawa couplings Yi^^ given in (14. 4p (solid 
line) and its two- heavy- neutrino mixing approximation (dashed line) where we set Rai3 = 0. 
Note that the two-heavy-neutrino mixing approximation can differ from the full calculation 
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Figure 8: The baryon asymmetry rj^ (blue contours) and the branching ratio B{t — t- 
/i7) (red contours) as functions of ki and k,2 the ReL model with = 120 GeV, 
considering a normal light neutrino mass spectrum. The following choice of the remaining 
parameters: 71 = Svr/S, 72 = 7r/2, 0i = vr, 02 = 0, Re{a) > (left panel); 71 = tt/2, 
72 = Svr/S, 01 = 0, 02 = TT, Re{a) > (right panel). The neutrino oscillation parame- 
ters are chosen at their best fit values but with sin^ 613 = 0.033 at its 2a upper limit. The 
blue shaded regions denote the parameter space where the baryon asymmetry is larger than 
rj'^^ = 6.2 X 10"^'^. The green shaded areas labeled as '6m.'^ = 0.25' and '0.10' indicate 
the parameter space where the inversion of the light-neutrino mass matrix is violated at the 
25% and 10% levels, respectively. 



by up to one order of magnitude. In addition, Figure E] shows the baryon asymmetry rjs 
calculated by omitting the RlS-subtracted collision terms in the BE f l4.27p (dotted line), 
or equivalently by taking the parameter Kr defined in f l4.36p equal to 1 in the BE f l4.33p . 
Such a simplification may reduce the predicted values for by as much as 60%. 

Figures [6] and [7] present numerical estimates, for a R/iL scenario with normal and 
inverted light neutrino mass spectra, respectively. We see that the baryon asymmetry rjs 
exhibits a similar dependence on /ti and K2 as in the RrL scenario. Due to the large e- 
and r- Yukawa couplings present in the R/iL model, the largest LEV rate can be observed 
in the r — e transitions, e.g. in the LEV process r — )■ 67, with B{t — )■ 67) ~ 10"^°. 
Since the current experimental sensitivity to this process is Bcxp{T — )■ 67) ~ 10"'^ which is 
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Figure 9: The same as in Figurel^ but for an inverted light neutrino mass spectrum and 
the following choice of parameters: 71 = Svr/S, 72 = vr/2, 0i = 0, 02 = ti"; Re{a) < (left 
panel); 71 = 7r/2, 72 = Svr/S, 0i = 0, 02 = 0, i?e(a) > (right panel). 

not expected to increase by more than one order of magnitude in the foreseeable future, 
it would be difficult to probe the parameter space of the R/xL scenario compatible with 
observable BAU. On the other hand, the — e transitions in the R/iL model, although 
being proportional to max(Ki, ^2)^0^ and so smaller than the predictions obtained in RrL 
scenario, are still sizeable enough to produce a /i — > e conversion rate of -8^6(22^1) ~ 
2 X 10"^'' (for Am^3 > 0) and 5 x 10~^^ (for Am^3 < 0). As the parameters a and b are 
approximately inversely proportional to ki,2, these values are largely independent of ki 2 and 
apply to the whole (ki, ^2) parameter plane, as depicted in Figures [6] and [3 Consequently, 
this scenario could be probed at a future /x — )■ e conversion experiment with a sensitivity 
of ~ IQ-^^ (COMET, mu2e), or ~ lO^^*^ (PRISM) [21j. 

Finally, Figures [8] and [9] display numerical estimates, for a ReL scenario with normal 
and inverted light neutrino mass spectra, respectively. Our results are quite analogous to 
the R/iL case. Correspondingly, the largest LFV rate is obtained for the process r — > fi'j. 
However, successful ReL requires B{t — )■ fi'-f) ~ 10~^^, which is far beyond the reach of the 
next generation experiments. In analogy to the KfiL case, the rates for coherent fi ^ e 
conversion in nuclei are found to be sizeable. Specifically, we obtain i?^e(22Ti) ^ 3 x 10~^^ 
(for Am^3 > 0) and 7 x 10~^^ (for Am^3 < 0). Interestingly enough, these rates are well 
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within reach of the proposed PRISM experiment [2T] . 



6 Conclusions 

We have analyzed minimal low-scale seesaw scenarios of resonant leptogenesis and studied 
their potential implications for observables of charged LFV, such as /i — ?■ 67 and /i — )■ e 
conversion in nuclei. We have considered three physically interesting flavour realisations 
of resonant leptogenesis, where the observed BAU originates from an individual r-, fi- 
or e-lepton-number asymmetry which gets resonantly enhanced via the out-of-equilibrium 
decays of nearly degenerate heavy Majorana neutrinos. 

By means of approximate lepton-fiavour symmetries, we have been able to construct 
viable and natural models of R£L compatible with universal right-handed neutrino masses 
at the GUT scale, where the required heavy-neutrino mass splittings are generated via 
RG effects. Particular attention has been paid that the effective resummation method 
introduced in |7j and used in our study to compute the resonantly enhanced lepton asym- 
metries respects the Nanopoulos- Weinberg no-go theorem [5l] in the L-conserving limit of 
the theory. Specifically, we have checked that the leptonic asymmetries Sai given in ( 14. Sp 
and (14.91) vanish in all parametrically possible L-conserving limits of the BIL scenarios. In 
agreement with earlier studies p^[H]. we find that at least three heavy Majorana neutrinos 
are required, in order to potentially have both successful leptogenesis and experimentally 
testable rates for LFV processes, such as /i — )■ 67 and /x — ?• e conversion in nuclei. 

We have found that the heavy Majorana neutrinos in R£L scenarios can be as light 
as 100 GeV, whilst their couplings to two of the charged leptons may be large so as to lead 
to LFV effects that could be tested by the MEG and the COMET/PRISM experiments. 
Specifically, in the RrL model with a normal light neutrino mass hierarchy, there is a 
sizeable model parameter space with successful leptogenesis and large LFV process rates, 
with B{fi — )■ 67) ~ 10^^^. This prediction is largely independent of sin^ 6'i3 and the other 
light neutrino oscillation parameters. On the other hand, in the RrL model with inversely 
hierarchical light neutrinos, — )■ 67) is linearly proportional to sin^6'i3, and can be 
enhanced by more than one order of magnitude compared to the normal hierarchy case 
for sin^ 613 close to its upper experimental 2cr limit. Unfortunately, the generated baryon 
asymmetry r/B is suppressed in this scenario, and to test the viable parameter space for 
successful leptogenesis would require an experiment for — )■ e conversion in nuclei which 
is sensitive to i?^e ~ 10~^''-10~^®. This feature is also quite generic for the case of the R/xL 
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and ReL models, where /i — ?■ e flavour transitions are suppressed by the small r- Yukawa 
couplings Ki 2- In all R£L models, charged LFV in the r-lepton sector turns out to be at 
least 6 orders of magnitude beyond the current experimental sensitivity, as the predicted 
branching ratios are B{t — )■ 67, /i7) < 10~^^ in parameter regions required for successful 
leptogenesis. 

Further studies will be needed to analyze the full range of theoretical, phenomeno- 
logical and cosmological implications of the three different universal models of ReL, R/xL 
and RrL. For instance, the consideration of thermal effects [58] may provide a significant 
improvement on the standard framework of classical BEs adopted in the present analysis. 
An equally significant issue is whether our minimal RL models can account for the well- 
known problem of cold Dark Matter (CDM) in the Universe. An obvious solution would be 
to consider supersymmetric versions of R£L scenarios and study the relic abundance 

of the lightest stable SUSY particle, which could be a thermal right-handed sneutrino [61] . 
Alternatively, one may consider scale-invariant extensions of the SM with right-handed 
neutrinos [621163] . which are minimally augmented with one complex singlet scalar field 
whose one-loop induced VEV can naturally explain the origin of the electroweak-scale 
mass of the heavy Majorana neutrinos. It has been shown recently [63] that a minimal 
Z4-symmetric variant of these models can stay perturbative up to the Planck scale, as 
well as provide a CDM candidate through the so-called Higgs-portal mechanism [Bl]. It is 
therefore rather motivating to perform a dedicated analysis of observing electroweak-scale 
heavy Majorana neutrinos within the specific context of the R£L scenarios studied here, 
through their possible like-sign dilepton [26|[65] and/or trilepton [66] signatures at the LHC 
or at other future high-energy colliders. 
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